projectUrl <- "https://github.com/GranderLab/miR34a_asRNA_project/raw/master/inst/"
dataUrl <- "https://github.com/GranderLab/miR34a_asRNA_project/raw/master/data"

Figure 1a

Introduction

An antisense transcript to the miR34a HG was detected using the UCSC genome browser.

Methods

The USCS genome browser was utilized for the bioinformatic evaluation of antisense transcription utilizing the Human EST and ENCODE RNAseq tracks.

Results

url <- "miR34aLocus/miR34aLocus_hg38.png"
knitr::include_graphics(file.path(projectUrl, url))


A schematic picture of the miR34a locus from UCSC genome browser (hg38) including miR34a HG and mature miR34a, and LINC01759. H3K4me3 ChIP-seq data and conservation over the locus is also shown.

Conclusion

GENCODE annotation supports the presence of a transcript orinionating in the antisense orientation from the miR34a gene.

Figure 1b

Introduction

We wished to examine expression in different tissue types and, therefore, utilized a panel of common cell lines to detect miR34a asRNA expression, as well as, miR34a HG. We included both p53+/+ and p53-/- cell lines in the cell line panel due to the fact that miR34a is a known downstream target of p53.

Methods

Cells were cultured with the appropriate growth medium and seeded at a desity between 50-75% confluency. The following day, RNA was extracted using the RNeasy mini kit (Qiagen) and subsequently treated with DNase (Ambion Turbo DNA-free, Life Technologies). 500ng RNA was used for cDNA synthesis using MuMLV (Life Technologies) and a 1:1 mix of oligo(dT) and random nanomers. PCR was run for 35 or 38 or 16 cycles for miR34a asRNA, miR34a HG, and beta-actin, respectivley, with 58˚C annealing temperature. PCR products were analyzed on a 2% agarose gel.

Primers
primers <- data.frame(
    name=c(
        "miR34a asRNA F1",
        "miR34a asRNA R1",
        "miR34aHG_F",
        "miR34aHG_R",
        "B-actin_F",
        "B-actin_R"
    ),
    sequence=c(
        "AGC GGC ATC TCC TCC ACC TGA AA",
        "TTG CCT CGT GAG TCC AAG GAG AAT",
        "TCT GCT CCA GTG GCT GAT GAG AAA",
        "GTT CAC TGG CCT CAA AGT TGG CAT",
        "AGG TCA TCA CCA TTG GCA ATG AG",
        "CTT TGC GGA TGT CCA CGT CA"
    )
)
primers
name sequence
miR34a asRNA F1 AGC GGC ATC TCC TCC ACC TGA AA
miR34a asRNA R1 TTG CCT CGT GAG TCC AAG GAG AAT
miR34aHG_F TCT GCT CCA GTG GCT GAT GAG AAA
miR34aHG_R GTT CAC TGG CCT CAA AGT TGG CAT
B-actin_F AGG TCA TCA CCA TTG GCA ATG AG
B-actin_R CTT TGC GGA TGT CCA CGT CA

Results

url <- "miR34aASExpressionCellPanel/cellPanelLarge.png"
knitr::include_graphics(file.path(projectUrl, url))


Semi-quantitative PCR data from the screening of a panel of cell lines. * Indicates that, although the cell line is p53wt, other mechanisms are present which inhibit p53 function.

Conclusions

We were able to detect co-expression of both miR34a HG and miR34a AS in a variety of cell lines, including HCT116 human colon cancer cells and 293T human embryonic kidney cells. A correlation was observed between both miR34a HG and miR34a asRNA expression and p53 status, with p53+/+ cell lines tending to have higher expression of both transcripts.

Figure 1c

Insert code here from Erik and write Intro Methods etc

Results

url <- "TCGAcorrelation/miR34a_RP3-510D11.2_diploid_only_edited.png"
knitr::include_graphics(file.path(projectUrl, url))


Figure 1d

Introduction

In order to define the 3’ transcription stop site for the miR34a asRNA, 3’-RACE was performed.

Methods

All cell lines were cultured at 5% CO2 and 37° C with U2OS cells grown in McCoy’s 5a (Life Technologies). All growth mediums were supplemented with 10% heat-inactivated FBS and 50 μg/ml of streptomycin and 50 μg/ml of penicillin. U2OS cell RNA was polyA-tailed using yeast polyA polymerase after which cDNA was synthesized using oligo(dT) primers. Nested-PCR was performed first using a forward primer in miR34a asRNA exon 1 (miR-34a asRNA F10 forward primer - ACGCGTCTCTCCAGCCCGGGAT) and a tailed oligo(dT) primer (poly T T7-FAM) followed by a second PCR using an alternate miR34a asRNA exon 1 primer (miR-34a asRNA F1 forward primer - AGCGGCATCTCCTCCACCTGAAA) and a reverse primer binding to the tail of the previously used oligo(dT) primer (FAM primer). PCR products were gel purified and cloned the Strata Clone Kit (Agilent Technologies), and sequenced.

Results

url <- "3primeRACE/results_hg38.png"
knitr::include_graphics(file.path(projectUrl, url))


3’-RACE sequencing results displayed in the UCSC genome browser together with the annotated miR34a asRNA transcript from GENCODE.

Conclusions

Sequencing of the resulting cloned cDNA indicated the transcripts 3’ transcription termination site to be 125bp upstream of the RP3-510D11.2 transcript’s annotated termination site.

Figure 1e

Introduction

We desired to determine the 5’-start site for the miR34a asRNA.

Methods

To accomplish this we performed the primer walk assay. All cell lines were cultured at 5% CO2 and 37° C with HEK293T cells cultured in DMEM high glucose (Hyclone). All growth mediums were supplemented with 10% heat-inactivated FBS and 50 μg/ml of streptomycin and 50 μg/ml of penicillin. RNA was extracted using the RNeasy mini kit (Qiagen) and subsequently treated with DNase (Ambion Turbo DNA-free, Life Technologies). 500ng RNA was used for cDNA synthesis using MuMLV (Life Technologies) and a 1:1 mix of oligo(dT) and random nanomers. PCR was performed for 35 cycles at 62°C annealing temperature and the resulting products were analyzed on a 2% agarose gel.

Primers
primers <- data.frame(
    name=c(
        "miR34aAS_F10",
        "miR34aAS_F11",
        "miR34aAS_F12",
        "miR34aAS_F13",
        "miR34aAS_F14",
        "miR34aAS_F15",
        "miR34aAS_R1"
    ),
    sequence=c(
        "ACG CGT CTC TCC AGC CCG GGA T",
        "ATC TGC GTG GTC ACC GAG AAG CA",
        "CGC ACG GAC TGA GAA ACA CAA G",
        "ACG GAG GCT ACA CAA TTG AAC AGG",
        "AGG GAA GAA AGA ACT AGC CGA GCA",
        "CAT TTG CTG CAA TAT CAC CGT GGC",
        "TTG CCT CGT GAG TCC AAG GAG AAT"
    )
)
primers
name sequence
miR34aAS_F10 ACG CGT CTC TCC AGC CCG GGA T
miR34aAS_F11 ATC TGC GTG GTC ACC GAG AAG CA
miR34aAS_F12 CGC ACG GAC TGA GAA ACA CAA G
miR34aAS_F13 ACG GAG GCT ACA CAA TTG AAC AGG
miR34aAS_F14 AGG GAA GAA AGA ACT AGC CGA GCA
miR34aAS_F15 CAT TTG CTG CAA TAT CAC CGT GGC
miR34aAS_R1 TTG CCT CGT GAG TCC AAG GAG AAT



Results

url <- "primerWalk/schematic_hg38.png"
knitr::include_graphics(file.path(projectUrl, url))


A schematic picture of the miR34a asRNA locus including the location of the primers used in the primer walk assay.

url <- "primerWalk/gel.png"
knitr::include_graphics(file.path(projectUrl, url))


The semi-quantitative PCR results from the primer walk assay performed using HEK293T cells.

Conclusion

It appears as though the miR34a asRNA transcript start site is somewhat upstream of the annotated start site. Although there is some product present for the F12 primer, the majority of transcription seems to be originating from the F11 primer. This indicates that the start site is between 71-224bp upstream of the annotated start site.

Figure 1f and supplementary figure ??????????????

Introduction

We utilized a bioinformatic approach to evaluate the coding potential of the miR34a asRNA transcript. The Coding-potential assessment tool uses a linear regression model to evaluate coding-potential by examining ORF length, ORF coverage, Fickett score and hexamer score. We further confirmed these results using the Coding-potential Calculator which utilizes a support based machine-based classifier and accesses an alternate set of discriminatory features.

Methods

Protein-coding capacity was evaluated using the Coding-potential Assessment Tool and Coding-potential Calculator with default settings. Transcript sequences for use with Coding-potential Assessment Tool were downloaded from the UCSC genome browser using the following IDs: HOTAIR (ENST00000455246.1), XIST (ENST00000429829.1), β-actin (ENST00000331789.5), Tubulin (ENST00000427480.1), and MYC (ENST00000377970). Transcript sequences for use with Coding-potential Calculator were downloaded from the UCSC genome browser using the following IDs: HOTAIR (uc031qho.1), β-actin (uc003soq.4).

Results

Coding Potential Assesment Tool
url <- "codingPotential/CodingPotentialAssessmentTool/results.txt"
cpat <-read.table(url(file.path(projectUrl, url)), header = TRUE, sep = "\t")
cpat
sequenceName RNAsize ORFsize FicketScore hexamerScore codingProbability codingLabel
miR34a asRNA 795 285 0.801 0.258 0.301 no
HOTAIR 918 144 0.919 0.268 0.126 no
XIST 19280 411 0.663 -0.128 0.027 no
β-actin 1917 1128 1.350 0.698 1.000 yes
Tubulin 2632 1119 1.277 0.515 1.000 yes
MYC 2345 1320 1.223 0.484 1.000 yes


Coding potential analysis results from the Coding-potential Assessment Tool including miR34a asRNA and two characterized non-coding transcripts (HOTAIR and XIST) and 3 known coding transcripts (β-actin, tubulin, and MYC).

Coding Potential Calculator
url <- "codingPotential/CodingPotentialCalculator/results.txt"
cpc <- read.table(url(file.path(projectUrl, url)), header = TRUE, sep = "\t")
cpc
Transcript.name Coding Coding.Potential.Score Hit.number Hit.score Frame.score Length Coverage…. Log.Odds.score Type
HOTAIR noncoding -1.187350 0 0.0000 0.00000 319 13.43 32.91 full
miR34a asRNA noncoding (weak) -0.883468 0 0.0000 0.00000 175 21.89 42.47 full
β-actin coding 13.662000 250 181.3849 34.26889 1167 62.82 229.54 full


Coding potential analysis results using the Coding-potential calculator for miR34a asRNA, HOTAIR, and β-actin.



Conclusions

Results indicated that miR34a asRNA has a similar lack of coding capacity to the known non-coding transcripts HOTAIR and XIST and differs greatly when examining these parameters to the known coding transcripts beta-actin, tubulin, and MYC. However, to fully evaluate coding potential methods such as mass spectrometry or ribosome profiling must be used.

Figure 2a

Introduction

miR34a is a known downstream target of p53 and has been previously shown to exhibit increased expression with cellular apoptotic signaling. We hypothesized that miR34a asRNA may be regulated in a similar fashion whereby transcription is stimulated by activation of p53. To test this we treated p53+/+ and p53-/- HCT116 cells with increasing concentrations of the DNA damaging agent doxorubicin and monitored miR34a asRNA expression.

Methods

HCT116 cells were cultured in DMEM high modified (Hyclone, GE healthcare) supplemented with supplemented with 2mM L-glutamine, 50ug/ml Penicillin-Streptomycin and 10% Fetal Calf Serum.
25 x 10^5 HCT116 cells were plated in 6 well plates. 24 hours later media was exchanged and doxorubicin was added to a final concentration of 100, 200 or 500 ng/ml. Cells were harvested for RNA extraction 24 hours later using trypsin. RNA was extracted using Nucleospin RNA kit (Machery-Nagel Ref. 740955) according to manufacturer‟s protocol and DNase treated using Ambion Turbo DNA-free according to manufacturer‟s protocol (Life Technologies Ref. AM1907). cDNA was synthesized using ~500 ng RNA with M-MLV (Life Technologies Ref 28025013) and a mixture of oligo(dT) with nanomers in accordance with the manufacturer’s protocol. qPCR quantification was carried out using the PowerUp SYBR Green Master Mix (Thermo Fisher Scientific, Ref. A25777) on the CFX96 Touch Real-Time PCR Detection System: 50°C for 2 min, 95°C for 2min, and 95°C for 1 sec followed by 60°C for 30 sec repeated for 40 cycles. The following primers were used:

Primers
primers <- data.frame(
    name=c(
        "ß-actin Fwd",
        "ß-actin Rev",
        "miR34a HG_F",
        "miR34a HG_R",
        "miR34a_asF1",
        "miR34a_asR1"
        
    ),
    sequence=c(
        "AGGTCATCACCATTGGCAATGAG",
        "CTTTCGGGATGTCCACGTCA",
        "TCTGCTCCAGTGGCTGATGAGAAA",
        "GTTCACTGGCCTCAAAGTTGGCAT",
        "AGCGGCATCTCCTCCACCTGAAA",
        "TTGCCTCGTGAGTCCAAGGAGAAT"
        
    )
)
primers
name sequence
ß-actin Fwd AGGTCATCACCATTGGCAATGAG
ß-actin Rev CTTTCGGGATGTCCACGTCA
miR34a HG_F TCTGCTCCAGTGGCTGATGAGAAA
miR34a HG_R GTTCACTGGCCTCAAAGTTGGCAT
miR34a_asF1 AGCGGCATCTCCTCCACCTGAAA
miR34a_asR1 TTGCCTCGTGAGTCCAAGGAGAAT


<
/br>

Results

Perform ddct analysis and statistical testing.

#download and format data
url <- "HCT116p53null.rda"
repmis::source_data(file.path(dataUrl, url))
## Downloading data from: https://github.com/GranderLab/miR34a_asRNA_project/raw/master/data/HCT116p53null.rda
## SHA-1 hash of the downloaded data file is:
## 2f3495af60d2b38454ecad092a9da798fafab11d
## [1] "."
data <- .

#calculate dct
groups <- c("Condition", "Treatment", "Biological Replicate", "Cell line", "gene")
gois <- c("miR34a asRNA", "miR34a HG")

data <- technicalMeans(data, groups) %>%
  dct(., gois, Actin, groups)

#calculate ddct
groups[groups == "gene"] <- "GOI"
logical <- tibble(
  Treatment = data$Treatment == 0
)

data <- data %>%
  ddct(., logical, groups) %>%
  folds(.)

#calculate stats
stats <- calcStats(data, Condition, "p53-wt", groups) %>%
  pFormat(.)

#setup plotting variables
stats <- stats %>%
  group_by(GOI, Treatment) %>%
  mutate(
    max = max(CI95h),
    x = case_when(
      Treatment ==   0 ~ 1 + 0.25,
      Treatment == 100 ~ 2 + 0.25,
      Treatment == 200 ~ 3 + 0.25,
      Treatment == 500 ~ 4 + 0.25
    ),
    xend = case_when(
      Treatment ==   0 ~ 1 - 0.25,
      Treatment == 100 ~ 2 - 0.25,
      Treatment == 200 ~ 3 - 0.25,
      Treatment == 500 ~ 4 - 0.25
    )
  ) %>%
  ungroup()

Plot the results.

#plot for Rmarkdown
ggplot(data = NULL) +
    geom_violin(
        data = data,
        aes(x = Treatment, y = log2fold, fill = Condition),
        trim = TRUE,
        alpha = 0.3,
        color = NA,
        scale = "width",
        show.legend = TRUE
    ) +
    geom_point(
        data = stats,
        aes(x = Treatment, y = mean, group = Condition), 
        position = position_dodge(width = 0.9), 
        size = 4
    ) +
    facet_grid(. ~ GOI) +
    geom_linerange(
        data = stats,
        aes(x = Treatment, ymin = CI95l, ymax = CI95h, group = Condition),
        colour = "black",
        position = position_dodge(width = 0.9),
        show.legend = FALSE
    ) +
    geom_label(
        data = stats,
        aes(x = Treatment, y = max + 0.22, label = pFormat),
        label.size = 0,
        label.padding = unit(0.01, "lines"),
        show.legend = FALSE,
        fill = "white",
        size = 5,
        family = "Arial Unicode MS",
        position = position_dodge(width = 0.9)
    ) +
    geom_segment(
        data = stats,
        aes(x = x, y = max + 0.1, xend = xend, yend = max + 0.1),
        colour="grey43",
        show.legend = FALSE
    ) +
    scale_y_continuous(
        "\u0394\u0394Ct",
        breaks=c(-1,0,1,2,3,4,5),
        sec.axis = sec_axis(
            ~2^., 
            name = expression(2^{- ~ Delta ~ Delta ~ "Ct"}), 
            breaks=c(2^-1,2^0,2^1,2^2,2^3,2^4,2^5)
        )
    )+
    labs(
      x="Treatment",
      title="miR34a asRNA regulation by p53 induction",
      caption=
            "Monitoring the effects of 24 hours doxorubicin treatment on miR34a asRNA and HG in P53+/+ and p53-/- HCT116 cells."
    )+
    theme_few()+
    scale_fill_ptol()+
    scale_colour_ptol()+
    #values=c("#014d64", "#6794a7")
    theme(
        plot.title = element_text(
            hjust = 0.5,
            face="bold",
            size=20
        ),
        plot.caption = element_text(
            hjust=0, 
            margin=margin(t=15),
            family="Arial Unicode MS"
        ),
        legend.position = "top",
        legend.title=element_text(size=17),
        legend.text=element_text(size=15),
        axis.title=element_text(size=17),
        axis.title.x=element_text(margin=margin(t=10)),
        axis.text=element_text(size=15),
        strip.text.x=element_text(size=17)
    )+
    guides(fill = guide_legend(title = "Cell line", override.aes = list(size=3)))



Conclusions

We observed increasing miR34a HG and asRNA expression with increasing doses of doxorubicin indicating that these two transcripts are co-regulated by p53. Although p53+/+ cells showed a 25-fold increase over untreated cells in miR34a AS expression at the highest tested doxorubicin concentration, this effect was largely abrogated, 11 fold over untreated, in p53-/- cells, indicating that p53 is a major regulator of miR34a asRNA expression.

Figure 2b

Introduction

miR34a is a known downstream target of p53 and has been previously shown to exhibit increased expression with cellular apoptotic signaling. We hypothesized that miR34a asRNA may be regulated in a similar fashion whereby transcription is stimulated by activation of p53. To test this we treated HCT116 and HEK293t cells with 200 ng/ml of the DNA damaging agent doxorubicin for 24 hours and monitored miR34a asRNA expression.

Methods

All cell lines were cultured at 5% CO2 and 37° C with HEK293T cells cultured in DMEM high glucose (Hyclone) and HCT116 cells in McCoy’s 5a (Life Technologies). All growth mediums were supplemented with 10% heat-inactivated FBS and 50 μg/ml of streptomycin and 50 μg/ml of penicillin. Cells were plated at 300,000 cells per well in a 6-well plate and cultured overnight. The following day cells were treated with 0, 100, 200, or 500 ng/ml doxorubicin for 24hrs. RNA was extracted using the RNeasy mini kit (Qiagen) and subsequently treated with DNase (Ambion Turbo DNA-free, Life Technologies). 500ng RNA was used for cDNA synthesis using MuMLV (Life Technologies) and a 1:1 mix of oligo(dT) and random nanomers. QPCR was carried out using KAPA 2G SYBRGreen (Kapa Biosystems) using the Applied Biosystems 7900HT machine with the cycling conditions: 95 °C for 3 min, 95 °C for 3 s, 60 °C for 30 s.

primers <- data.frame(
    name=c(
        "ß-actin Fwd",
        "ß-actin Rev",
        "miR34a HG_F",
        "miR34a HG_R",
        "miR34a_asF1",
        "miR34a_asR1"
        
    ),
    sequence=c(
        "AGGTCATCACCATTGGCAATGAG",
        "CTTTCGGGATGTCCACGTCA",
        "TCTGCTCCAGTGGCTGATGAGAAA",
        "GTTCACTGGCCTCAAAGTTGGCAT",
        "AGCGGCATCTCCTCCACCTGAAA",
        "TTGCCTCGTGAGTCCAAGGAGAAT"
        
    )
)
primers
name sequence
ß-actin Fwd AGGTCATCACCATTGGCAATGAG
ß-actin Rev CTTTCGGGATGTCCACGTCA
miR34a HG_F TCTGCTCCAGTGGCTGATGAGAAA
miR34a HG_R GTTCACTGGCCTCAAAGTTGGCAT
miR34a_asF1 AGCGGCATCTCCTCCACCTGAAA
miR34a_asR1 TTGCCTCGTGAGTCCAAGGAGAAT

Results

url <- "HctHekDox.rda"
repmis::source_data(file.path(dataUrl, url))
## Downloading data from: https://github.com/GranderLab/miR34a_asRNA_project/raw/master/data/HctHekDox.rda
## SHA-1 hash of the downloaded data file is:
## 54d24343591ee15b1c076723db38cb02b5142b4e
## [1] "."
data <- .

#calculate dct
groups <- c("Treatment", "Biological Replicate", "Cell line", "gene")
gois <- c("miR34a asRNA", "miR34a HG")

data <- technicalMeans(data, groups) %>%
  dct(., gois, Actin, groups)

#calculate ddct
groups[groups == "gene"] <- "GOI"
logical <- tibble(
  `Biological Replicate` = pull(data, `Biological Replicate`) == 1,
  Treatment = pull(data, Treatment) == "untreated"
)

data <- data %>%
  ddct(., logical, groups) %>%
  folds(.)

#calculate stats
stats <- calcStats(data, Treatment, "untreated", groups) %>%
  pFormat(.)
#plot for Rmarkdown
ggplot(data = NULL) +
  geom_violin(
    data = data,
    aes(x = Treatment, y = log2fold, fill = GOI),
    trim = TRUE,
    alpha = 0.3,
    color = NA,
    scale = "width",
    show.legend = FALSE
  ) +
  geom_point(
    data = stats,
    aes(x = Treatment, y = mean, group = GOI),
    size = 4,
    position = position_dodge(width = 0.9)
  ) +
  facet_grid(. ~ `Cell line`) +
  geom_linerange(
    data = stats,
    aes(x = Treatment, ymin = CI95l, ymax = CI95h, group = GOI),
    colour = "black",
    position = position_dodge(width = 0.9),
    show.legend = FALSE
  ) +
  geom_label(
    data = stats,
    aes(x = Treatment, y = CI95h + 0.15, label = pFormat, group = GOI),
    label.size = 0,
    label.padding = unit(0.01, "lines"),
    show.legend = FALSE,
    fill = "white",
    size = 5,
    family = "Arial Unicode MS",
    position = position_dodge(width = 0.9)
  )+
  scale_y_continuous(
    "\u0394\u0394Ct",
    breaks=c(-1,0,1,2,3,4,5)
  )+
  labs(
    x = "Treatment",
    title = "miR34a asRNA response to p53 induction",
    caption = 
            "Q-PCR showing miR34a asRNA and miR34a HG levels in HCT116 and HEK293t cells after treatment with 200 ng/ml doxorubicin for 24hrs.\nError bars (black vertical) indicate 95% CI and means are represented \nby black dots."
  )+
  theme_few()+
  scale_fill_ptol()+
  scale_colour_ptol()+
  theme(
    plot.title = element_text(
      hjust = 0.5,
      face = "bold",
      size = 20
    ),
    plot.caption = element_text(
      hjust = 0, 
      margin = margin(t = 15),
      family = "Arial Unicode MS"
    ),
    legend.position = "top",
    legend.title = element_text(size = 17),
    legend.text = element_text(size = 15),
    axis.title = element_text(size = 17),
    axis.title.x = element_text(margin = margin(t = 10)),
    axis.text = element_text(size = 15),
    strip.text.x = element_text(size = 17)
  )+
  guides(color = guide_legend(title = "Gene", override.aes = list(size = 3)))

Conclusions

We observed increases in miR34a HG and asRNA expression upon doxorubicin treatment indicating that these two transcripts are co-regulated.

Figure 2c

Introduction

It is likely, due to the head-to head orientation of miR34a HG and asRNA, that transcription initiation may be activated from a single promoter in a bi-directional manner. To investigate whether miR34a HG and asRNA are transcribed from the same promoter as convergent transcripts, we cloned the miR34a HG promoter, including the p53 binding site, into a luciferase/renilla dual reporter vector which we hereafter refer to as p1. The p1 sequence was previously published in Raver-Shapira, N., et al., Transcriptional activation of miR-34a contributes to p53-mediated apoptosis. Mol Cell, 2007. 26(5): p. 731-43. Due to the observance of additional conservation sites downstream of the previously published miR34a HG promoter, we also cloned a longer “promoter region”, extending 251 additional base pairs into the miR34a HG into the same bi-directional reporter vector, hereafter refered to as p2.

Methods

All cell lines were cultured at 5% CO2 and 37° C with HEK293T cells cultured in DMEM high glucose (Hyclone) and HCT116 cells in McCoy’s 5a (Life Technologies). All growth mediums were supplemented with 10% heat-inactivated FBS and 50 μg/ml of streptomycin and 50 μg/ml of penicillin. Cells were plated at 10,000 cells per well in a 96-well plate with a white bottom and cultured overnight. The following day cells were co-transfected with 10ng of empty, p1, or p2 plasmid and GFP using the standard lipofektamine 2000 (Life Technologies) protocol. The expression of GFP and luminescence was measured 24 h post transfection by using the Dual-Glo Luciferase Assay System (Promega) and detected by the GloMax-Multi+ Detection System (Promega). The expression of luminescence was normalized to GFP.

Constructs
library(printr)
ps <- data.frame(
    name=c("p1", "p2"),
    sequence=c("GCG CCC TGC CTG GCC CCC ACC TGG TCC TCT TTC CTT TTC AGG TGG AGG AGA TGC CGC TGT
              CCC GTC GGT CTG GGG ACA GCC CAG CTC CCC GGA TCC CGG GCT GGA GAG ACG CGT CGC GGC   
              CCC GGG GCC TGG TGG CAC GAG CAG GAA GGA GGA CCC GGC GGC GGG CTC TGC CTG GGC TTG 
              CCT GGG CTT GTT CCG AGC CGG GCT GCT TCT CGG TGA CCA CGC AGA TCG GGG GCA TTT GGA
              GAT TTT GCG GGA GTC CTG CAG CCA AGC TCC GGG GCA GGA GAG GCC TGG AAG CCT GCA CTA
              CCT GCT C",
              "GCG CCC TGC CTG GCC CCC ACC TGG TCC TCT TTC CTT TTC AGG TGG AGG AGA TGC CGC TGT 
              CCC GTC GGT CTG GGG ACA GCC CAG CTC CCC GGA TCC CGG GCT GGA GAG ACG CGT CGC GGC
              CCC GGG GCC TGG TGG CAC GAG CAG GAA GGA GGA CCC GGC GGC GGG CTC TGC CTG GGC TTG 
              CCT GGG CTT GTT CCG AGC CGG GCT GCT TCT CGG TGA CCA CGC AGA TCG GGG GCA TTT GGA 
              GAT TTT GCG GGA GTC CTG CAG CCA AGC TCC GGG GCA GGA GAG GCC TGG AAG CCT GCA CTA 
              CCT GCT CGC CCC GTC CCA GCA TGC ACC CAG GTA AAC GCT TGT GTT TCT CAG TCC GTG CGA 
              AAG TTT GCA AAG AAG GAG GCG GGA ACT AGA CCA ACA ACT TTA ATA ATC ATA GTC ATA ATA 
              ATG AAA ACC CTG TTC AAT TGT GTA GCC TCC GTA AGG GGA AGA ATT TCC AAG AAG CAA CAA
              AAT GTT CTT TGA TTT CAT AAT CTC TGG GGA GAA AGG ATT CCT TTT TTG GGG TTT TGC TCG
              GCT AGT TCT TTC TTC CCT")
)
ps
name sequence
p1 GCG CCC TGC CTG GCC CCC ACC TGG TCC TCT TTC CTT TTC AGG TGG AGG AGA TGC CGC TGT
CCC GTC GGT CTG GGG ACA GCC CAG CTC CCC GGA TCC CGG GCT GGA GAG ACG CGT CGC GGC
CCC GGG GCC TGG TGG CAC GAG CAG GAA GGA GGA CCC GGC GGC GGG CTC TGC CTG GGC TTG
CCT GGG CTT GTT CCG AGC CGG GCT GCT TCT CGG TGA CCA CGC AGA TCG GGG GCA TTT GGA
GAT TTT GCG GGA GTC CTG CAG CCA AGC TCC GGG GCA GGA GAG GCC TGG AAG CCT GCA CTA
CCT GCT C
p2 GCG CCC TGC CTG GCC CCC ACC TGG TCC TCT TTC CTT TTC AGG TGG AGG AGA TGC CGC TGT
CCC GTC GGT CTG GGG ACA GCC CAG CTC CCC GGA TCC CGG GCT GGA GAG ACG CGT CGC GGC
CCC GGG GCC TGG TGG CAC GAG CAG GAA GGA GGA CCC GGC GGC GGG CTC TGC CTG GGC TTG
CCT GGG CTT GTT CCG AGC CGG GCT GCT TCT CGG TGA CCA CGC AGA TCG GGG GCA TTT GGA
GAT TTT GCG GGA GTC CTG CAG CCA AGC TCC GGG GCA GGA GAG GCC TGG AAG CCT GCA CTA
CCT GCT CGC CCC GTC CCA GCA TGC ACC CAG GTA AAC GCT TGT GTT TCT CAG TCC GTG CGA
AAG TTT GCA AAG AAG GAG GCG GGA ACT AGA CCA ACA ACT TTA ATA ATC ATA GTC ATA ATA
ATG AAA ACC CTG TTC AAT TGT GTA GCC TCC GTA AGG GGA AGA ATT TCC AAG AAG CAA CAA
AAT GTT CTT TGA TTT CAT AAT CTC TGG GGA GAA AGG ATT CCT TTT TTG GGG TTT TGC TCG
GCT AGT TCT TTC TTC CCT
Primers
primers <- data.frame(
    name=c(
        "Luc_set_II_F",
        "Luc_set_II_R",
        "Renilla_pBiDir_F1",
        "Renilla_pBiDir_R1",
        "B-actin_F",
        "B-actin_R"
    ),
    sequence=c(
        "AAG ATT CAA AGT GCG CTG CTG",
        "TTG CCT GAT ACC TGG CAG ATG",
        "TAA CGC GGC CTC TTC TTA TTT",
        "GAT TTG CCT GAT TTG CCC ATA",
        "AGG TCA TCA CCA TTG GCA ATG AG",
        "CTT TGC GGA TGT CCA CGT CA"
    )
)
primers
name sequence
Luc_set_II_F AAG ATT CAA AGT GCG CTG CTG
Luc_set_II_R TTG CCT GAT ACC TGG CAG ATG
Renilla_pBiDir_F1 TAA CGC GGC CTC TTC TTA TTT
Renilla_pBiDir_R1 GAT TTG CCT GAT TTG CCC ATA
B-actin_F AGG TCA TCA CCA TTG GCA ATG AG
B-actin_R CTT TGC GGA TGT CCA CGT CA

Results

url <- "P1-HCTandHEK/locusSimple.png"
knitr::include_graphics(file.path(projectUrl, url))



url <- "P1-HCTandHEK/constructSimple.png"
knitr::include_graphics(file.path(projectUrl, url))



url <- "P1-HCTandHEK.rda"
repmis::source_data(file.path(dataUrl, url))
## Downloading data from: https://github.com/GranderLab/miR34a_asRNA_project/raw/master/data/P1-HCTandHEK.rda
## SHA-1 hash of the downloaded data file is:
## cec9138dee1f4183d9dbfc1f039a11b3024232fd
## [1] "."
data <- .

.calcP <- function(`Cell line`, gene, value, data) {
  bool1 <- pull(data, `Cell line`) == `Cell line`
  bool2 <- pull(data, gene) == gene
  bool3 <- pull(data, construct) == "Empty"
  bool <- bool1 & bool2 & bool3
  t.test(value, pull(data, value)[bool])$p.value
}

stats <- data %>%
  group_by(`construct`, `Cell line`, gene) %>%
  summarize(
    n = n(),
    mean = mean(value),
    CI95l = t.test(value)$conf.int[1],
    CI95h = t.test(value)$conf.int[2],
    pValue = .calcP(`Cell line`, gene, value, data)
  ) %>%
  ungroup() %>%
  pFormat(.)
data <- filter(data, construct %in% c("p1", "Empty"))
stats <- filter(stats, construct %in% c("p1", "Empty"))

#plot for Rmarkdown
ggplot(data = NULL) +
  geom_violin(
    data = data,
      aes(x = construct, y = value, fill = gene),
      trim = TRUE,
      alpha = 0.3,
      color = NA,
      scale = "width",
      show.legend = FALSE
  ) +
  geom_point(
    data = data,
    aes(x = construct, y = value, colour = gene),
    size = 4,
    position = position_dodge(width = 0.9),
    stroke = 0
  ) +
  facet_grid(. ~ `Cell line`) +
  geom_linerange(
    data = stats,
    aes(x = construct, ymin = CI95l, ymax = CI95h, group = gene),
    colour = "black",
    position = position_dodge(width = 0.9),
    show.legend = FALSE
  ) +
  geom_point(
    data = stats,
    aes(x = construct, y = mean, group = gene),
    colour = "black",
    position = position_dodge(width = 0.9),
    show.legend = FALSE
  ) +
  geom_label(
    data = stats,
    aes(x = construct, y = CI95h + 1900, label = pFormat, group = gene),
    label.size = 0,
    label.padding = unit(0.01, "lines"),
    show.legend = FALSE,
    fill = "white",
    size = 5,
    family = "Arial Unicode MS",
    position = position_dodge(width = 0.9)
  ) +
  labs(
    x = "Construct",
    y = "Activity (a.u.)",
    title = "miR34a asRNA regulation by miR34a promoter",
    caption = 
            "miR34a asRNA and miR34a HG levels in HCT116 and HEK293t cells after transfection with P1 and P2 constructs.\nLuciferase represents miR34a HG and renilla represents miR34a asRNA.\nError bars (black vertical) indicate 95% CI and means are represented \nby black dots."
    ) +
  theme_few() +
  scale_fill_ptol() +
  scale_colour_ptol() +
  theme(
    plot.title = element_text(
      hjust = 0.5,
      face = "bold",
      size = 20
    ),
    plot.caption = element_text(
      hjust = 0, 
      margin = margin(t = 15),
      family = "Arial Unicode MS"
    ),
    legend.position = "top",
    legend.title = element_text(size = 17),
    legend.text = element_text(size = 15),
    axis.title = element_text(size = 17),
    axis.title.x = element_text(margin = margin(t = 10)),
    axis.text = element_text(size = 15),
    strip.text.x = element_text(size = 17)
  ) +
  guides(color = guide_legend(title = "Gene", override.aes = list(size = 3)))

Conclusion

Upon transfection of this construct in two separate cell lines we observed increases in both luciferase and renilla indicating that miR34a HG and asRNA expression is regulated by a single promoter contained within the p1 construct. Transfection of the P2 construct showed to cause a reduction in luciferase (miR34a HG) levels but not renilla (miR34a AS) levels. This could be due to the presence of a STAT3 binding site in p2 which has been previously reported to negativley regulate miR34a HG.

Figure 2d

Introduction

Our previous results showed that miR34a HG and asRNA are simultaneously up-regulated in response to doxorubicin treatment leading us to consider the possibility that miR34a asRNA may be necessary for an appropriate miR34a response to apoptotic signaling. To explore this possibility we again utilized the p1 construct, co-transfecting it with an shRNA targeting renilla, and treating HCT116 cells with increasing levels of doxorubicin.

Methods

We utilized the P1 construct where the overlapping region of miR34a HG and miR34a AS is cloned with luciferase downstream of miR34a HG and renilla downstream of miR34a AS. The p1 sequence was previously published in Raver-Shapira, N., et al., Transcriptional activation of miR-34a contributes to p53-mediated apoptosis. Mol Cell, 2007. 26(5): p. 731-43. All cell lines were cultured at 5% CO2 and 37° C with HCT116 cells cultured in McCoy’s 5a (Life Technologies). All growth mediums were supplemented with 10% heat-inactivated FBS and 50 μg/ml of streptomycin and 50 μg/ml of penicillin. 2.5x10^5 HCT116wt cells were seeded in a 12-well plate. After 24hrs these were co-transfected with the P1 construct (25ng) and shRenilla2.1 (250ng) using lipo2000 (Life Technologies) standard protocol. 24hrs post-transfection, doxorubicin treatment was initiated using 0, 300, or 500ng/ml. 24hrs post-treatment, RNA was extracted using the RNeasy mini kit (Qiagen) and subsequently treated with DNase (Ambion Turbo DNA-free, Life Technologies). 500ng RNA was used for cDNA synthesis using MuMLV (Life Technologies) and a 1:1 mix of oligo(dT) and random nanomers. QPCR was carried out using KAPA 2G SYBRGreen (Kapa Biosystems) using the Applied Biosystems 7900HT machine with the cycling conditions: 95 °C for 3 min, 95 °C for 3 s, 60 °C for 30 s.

Constructs
library(printr)
p1 <- data.frame(
    name="p1",
    sequence="GCG CCC TGC CTG GCC CCC ACC TGG TCC TCT TTC CTT TTC AGG TGG AGG AGA TGC CGC TGT
              CCC GTC GGT CTG GGG ACA GCC CAG CTC CCC GGA TCC CGG GCT GGA GAG ACG CGT CGC GGC   
              CCC GGG GCC TGG TGG CAC GAG CAG GAA GGA GGA CCC GGC GGC GGG CTC TGC CTG GGC TTG 
              CCT GGG CTT GTT CCG AGC CGG GCT GCT TCT CGG TGA CCA CGC AGA TCG GGG GCA TTT GGA
              GAT TTT GCG GGA GTC CTG CAG CCA AGC TCC GGG GCA GGA GAG GCC TGG AAG CCT GCA CTA
              CCT GCT C"
)
p1
name sequence
p1 GCG CCC TGC CTG GCC CCC ACC TGG TCC TCT TTC CTT TTC AGG TGG AGG AGA TGC CGC TGT
CCC GTC GGT CTG GGG ACA GCC CAG CTC CCC GGA TCC CGG GCT GGA GAG ACG CGT CGC GGC
CCC GGG GCC TGG TGG CAC GAG CAG GAA GGA GGA CCC GGC GGC GGG CTC TGC CTG GGC TTG
CCT GGG CTT GTT CCG AGC CGG GCT GCT TCT CGG TGA CCA CGC AGA TCG GGG GCA TTT GGA
GAT TTT GCG GGA GTC CTG CAG CCA AGC TCC GGG GCA GGA GAG GCC TGG AAG CCT GCA CTA
CCT GCT C
shRNAs
shRNAs <- data.frame(
    name=c(
        "shRenilla2.1"
    ),
    sequence=c(
        "TAA CGG GAT TTC ACG AGG C"
    )
)
shRNAs
name sequence
shRenilla2.1 TAA CGG GAT TTC ACG AGG C
Primers
primers <- data.frame(
    name=c(
        "Luc_set_II_F",
        "Luc_set_II_R",
        "Renilla_pBiDir_F1",
        "Renilla_pBiDir_R1",
        "B-actin_F",
        "B-actin_R"
    ),
    sequence=c(
        "AAG ATT CAA AGT GCG CTG CTG",
        "TTG CCT GAT ACC TGG CAG ATG",
        "TAA CGC GGC CTC TTC TTA TTT",
        "GAT TTG CCT GAT TTG CCC ATA",
        "AGG TCA TCA CCA TTG GCA ATG AG",
        "CTT TGC GGA TGT CCA CGT CA"
    )
)
primers
name sequence
Luc_set_II_F AAG ATT CAA AGT GCG CTG CTG
Luc_set_II_R TTG CCT GAT ACC TGG CAG ATG
Renilla_pBiDir_F1 TAA CGC GGC CTC TTC TTA TTT
Renilla_pBiDir_R1 GAT TTG CCT GAT TTG CCC ATA
B-actin_F AGG TCA TCA CCA TTG GCA ATG AG
B-actin_R CTT TGC GGA TGT CCA CGT CA

Results

url <- "p1shRNAdox.rda"
repmis::source_data(file.path(dataUrl, url))
## Downloading data from: https://github.com/GranderLab/miR34a_asRNA_project/raw/master/data/p1shRNAdox.rda
## SHA-1 hash of the downloaded data file is:
## c41f628771e1915efb384c6483bc6cf6461d572e
## [1] "."
data <- .

#Luciferase = miR34a HG
#Renilla = miR34a AS

#calculate dct
groups <- c("treatment", "Biological Replicate", "shRNA", "gene")
gois <- c("Luciferase", "Renilla")

data <- data %>%
  technicalMeans(., groups) %>%
  dct(., gois, Actin, groups)

#calculate ddct
groups[groups == "gene"] <- "GOI"
logical <- tibble(
  `shRNA` = pull(data, shRNA) == "shRNA Control"
)

data <- data %>%
  ddct(., logical, groups) %>%
  folds(.)

#calculate stats
stats <- calcStats(data, shRNA, "shRNA Control", groups) %>%
  pFormat(.)

#prepare for plotting
data <- data %>%
  mutate(plotGOI = case_when(
    GOI == "Luciferase" ~ "  Luciferase\n(miR34a HG)",
    GOI == "Renilla" ~ "       Renilla\n(miR34a asRNA)"
  ))

stats <- stats %>%
  mutate(plotGOI = case_when(
    GOI == "Luciferase" ~ "  Luciferase\n(miR34a HG)",
    GOI == "Renilla" ~ "       Renilla\n(miR34a asRNA)"
  ))
#plot for rmarkdown
ggplot(data = NULL) +
  geom_violin(
    data = data,
    aes(x = treatment, y = log2fold, fill = plotGOI),
    trim=TRUE,
    scale="width",
    alpha=0.3,
    colour=NA,
    show.legend = FALSE
  ) +
  geom_point(
    data = data,
    aes(x = treatment, y = log2fold, colour = plotGOI),
    size = 4,
    position = position_dodge(width = 0.9)
  ) +
  facet_grid(. ~ shRNA, scale = "free") +
  geom_linerange(
    data = stats,
    aes(x = treatment, ymin = CI95l, ymax = CI95h, group = plotGOI),
    colour = "black",
    position = position_dodge(width = 0.9),
    show.legend = FALSE
  ) +
  geom_point(
    data = stats,
    aes(x = treatment, y = mean, group = plotGOI),
    colour = "black",
    position = position_dodge(width = 0.9),
    show.legend = FALSE
  ) +
  geom_label(
    data = stats,
    aes(x = treatment, y = CI95h + 0.1, label = pFormat, group = plotGOI),
    label.size = 0,
    label.padding = unit(0.01, "lines"),
    show.legend = FALSE,
    position = position_dodge(width = 0.9),
    fill = "white",
    size = 5,
    family = "Arial Unicode MS"
  ) +
  scale_y_continuous("\u0394\u0394Ct") +
  labs(
    x = "Doxorubicin (ng/ml)",
    title = "HCT116",
    caption =
            "miR34a (luficerase) and miR34a asRNA (renilla) in HCT116 cells which were co-transfected with the P1 construct and the specified shRNA after which they\nwere treated with the indicated dose of doxorubicin. Error bars (black vertical) indicate 95% CI and black points represent the mean."
  ) +
  theme_few() +
  scale_fill_ptol() +
  scale_color_ptol() +
  theme(
    plot.title = element_text(
      hjust = 0.5,
      face = "bold",
      size = 20
    ),
    plot.caption = element_text(
      hjust = 0, 
      margin = margin(t = 15),
      family = "Arial Unicode MS"
    ),
    legend.position = "top",
    legend.title = element_text(size = 17),
    legend.text = element_text(size = 15, hjust = 0.5),
    axis.title = element_text(size = 17),
    axis.title.x = element_text(margin = margin(t = 10)),
    axis.text = element_text(size = 15)
  ) + 
  guides(
    color = guide_legend(
      title = "Gene", 
      override.aes = list(size = 4)
    )
  )

Conclusion

Levels of luciferase were abrogated in the shRenilla-transfected cells compared to control shRNA-transfected cells. Collectively, these results indicate that miR34a asRNA positively regulates levels of miR34a HG and is crucial for a proper miR34a response to apoptotic stimuli.

Figure 3a

Introduction

We hypothesized that miR34a asRNA may function to mediate the levels of its sense gene, miR34a. To this end, we stably over-expressed miR34a asRNA, using lentiviral transduction, in 3 different p53-null cell lines and monitored the effects on miR34a.

Methods

miR34a asRNA cloning

RNA was harvested from U2OS cells using the RNeasy mini kit (Qiagen) and subsequently treated with DNase (Ambion Turbo DNA-free, Life Technologies). 500ng RNA was used for cDNA synthesis using MuMLV (Life Technologies) and a 1:1 mix of oligo(dT) and random nanomers. PCR was run using the miR34aAS_cloning_F4 and miR34aAS_cloning_Ex3_1 primers for 40 cycles at annealing temperature 58˚C and run on an agarose gel. Bands were excised and extracted using the QIAquick gel extraction kit (Qiagen). PCR products were cloned using the Strata Clone kit (Aligent) and sequenced (results below). The insert for the lentiviral construct (p-HIV-I2AG-U6) was subsequently sub-cloned from the sequenced plasmid using primers with PacI and NheI cleavage sites. The cloned transcript largley mirrors the GENCODE annotated transcript.

Cloning primers
library(printr)
primers <- data.frame(
    name=c(
        "miR34aAS_cloning_F4",
        "miR34aAS_cloning_Ex3_1"
    ),
    sequence=c(
        "ACGCGTCTCTCCAGCCCGGGAT",
        "AATGATGGCCGCAACTAATGACGG"
    )
)
primers
name sequence
miR34aAS_cloning_F4 ACGCGTCTCTCCAGCCCGGGAT
miR34aAS_cloning_Ex3_1 AATGATGGCCGCAACTAATGACGG
Sequencing results
construct <- data.frame(
    name="miR34a_asRNA",
    sequence="ACG CGT CTC TCC AGC CCG GGA TCC GGG GAG CTG GGC TGT CCC CAG ACC GAC GGG ACA GCG GCA
                TCT CCT CCA CCT GAA AAG GAA AGA GGA CCA GTT TGC AGG ACT CCG AAC TGG GCC CGC GAG ATC 
                TCC ACC TGC GCA AAA CGA AAG GGC GGA TTC TCC TTG GAC TCA CGA GGC AAC CGC TCC CCG GGG 
                TGA GAA CGG GGG ACT CAT TCC TCC GGC ACT GGG AGA AGA CGA TTC TTC AGG AGG AGG ACA GGG 
                AAG CGA ATG CTA CCC AGA TGT CTC AGT ATA CTG GCT CGC GGC ACA TCG GGC AAA TGA ACC TAT 
                CAG ATA ACA ACG GCA GAT CAG ATG CCT GAG CAT TCA GAA GCA ACA GCT GTG GAG CCC CCG TGG 
                GTT CAG AAG GCC TGG TTC CCG TCT CCA GAA GCC TGG CTC TCC TCC CTC CTG GGC CCA CTA CTT 
                TGG CTT CTT GTT CCT ACG TAC AAG GAG TTG CGA AGA AGG CAA CTC TTC CCC TCC CTG AAG CCA
                AAG GAA TGA AAC AGA CTA GGG CGG GAG AGG TGG CCA TCC GTC ATT AGT TGC GGC CAT CA"
)
construct
name sequence
miR34a_asRNA ACG CGT CTC TCC AGC CCG GGA TCC GGG GAG CTG GGC TGT CCC CAG ACC GAC GGG ACA GCG GCA
TCT CCT CCA CCT GAA AAG GAA AGA GGA CCA GTT TGC AGG ACT CCG AAC TGG GCC CGC GAG ATC
TCC ACC TGC GCA AAA CGA AAG GGC GGA TTC TCC TTG GAC TCA CGA GGC AAC CGC TCC CCG GGG
TGA GAA CGG GGG ACT CAT TCC TCC GGC ACT GGG AGA AGA CGA TTC TTC AGG AGG AGG ACA GGG
AAG CGA ATG CTA CCC AGA TGT CTC AGT ATA CTG GCT CGC GGC ACA TCG GGC AAA TGA ACC TAT
CAG ATA ACA ACG GCA GAT CAG ATG CCT GAG CAT TCA GAA GCA ACA GCT GTG GAG CCC CCG TGG
GTT CAG AAG GCC TGG TTC CCG TCT CCA GAA GCC TGG CTC TCC TCC CTC CTG GGC CCA CTA CTT
TGG CTT CTT GTT CCT ACG TAC AAG GAG TTG CGA AGA AGG CAA CTC TTC CCC TCC CTG AAG CCA
AAG GAA TGA AAC AGA CTA GGG CGG GAG AGG TGG CCA TCC GTC ATT AGT TGC GGC CAT CA
Lentivirus production and infection

HEK293T cells were transfected with viral and expression constructs using Lipofectamine 2000 (Life Technologies), after which viral supernatants were harvested 48 and 72 hours post-transfection. Viral particles were concentrated using PEG-IT solution (Systems Biosciences) according to the manufacturer’s recommendations. HEK293T cells were used for virus titration and GFP expression was evaluated 72hrs post-infection via flow cytometry after which TU/ml was calculated. Cell lines were subsequently infected with the lentiviral particals using a multuplicity of infection = 1. After 48 hours, cells were placed under mycophenolic acid selection and the selection process was monitored using GFP expression via flow cytometry. When GFP expression was present in >90% of the cells, cells were frozen for downstream analysis.

QPCR

RNA was extracted using the RNeasy mini kit (Qiagen) and subsequently treated with DNase (Ambion Turbo DNA-free, Life Technologies). 500ng RNA was used for cDNA synthesis using MuMLV (Life Technologies) and a 1:1 mix of oligo(dT) and random nanomers. QPCR was carried out using KAPA 2G SYBRGreen (Kapa Biosystems) using the Applied Biosystems 7900HT machine with the cycling conditions: 95 °C for 3 min, 95 °C for 3 s, 60 °C for 30 s. QPCR for miRNA expression analysis was performed according to the protocol for the TaqMan microRNA Assay kit (Life Technologies) with primer/probe sets purchased from the same company (TaqMan® MicroRNA Assay, hsa-miR-34a, human and Control miRNA Assay, RNU48, human) and the same cycling scheme as above.

QPCR primers
###Q-PCR primers
primers <- data.frame(
    name=c(
        "ß-actin Fwd",
        "ß-actin Rev",
        "miR34a_asF1",
        "miR34a_asR1"
        
    ),
    sequence=c(
        "AGGTCATCACCATTGGCAATGAG",
        "CTTTCGGGATGTCCACGTCA",
        "AGCGGCATCTCCTCCACCTGAAA",
        "TTGCCTCGTGAGTCCAAGGAGAAT"
        
    )
)
primers
name sequence
ß-actin Fwd AGGTCATCACCATTGGCAATGAG
ß-actin Rev CTTTCGGGATGTCCACGTCA
miR34a_asF1 AGCGGCATCTCCTCCACCTGAAA
miR34a_asR1 TTGCCTCGTGAGTCCAAGGAGAAT

Results

url <- "stableLineExpression.rda"
repmis::source_data(file.path(dataUrl, url))
## [1] "."
data <- .

#remove Skov exp 1 due to variability in technical replicates
data <- filter(data, `Biological Replicate` != 1 | `Cell line` != "Skov3")

#Have to divide the analysis due to the presence of several house 
#keeping genes.
gois1 <- c("miR34a asRNA", "Actin")
gois2 <- c("miR34a", "RNU48")
data1 <- filter(data, gene %in% gois1)
data2 <- filter(data, gene %in% gois2)

#calculate dct
groups <- c("Genetic mod", "Biological Replicate", "Cell line", "gene")

data1 <- data1 %>%
  technicalMeans(., groups) %>%
  dct(., gois1, Actin, groups) %>%
  mutate(GOI = as.character(GOI))

data2 <- data2 %>%
  technicalMeans(., groups) %>%
  dct(., gois2, RNU48, groups) %>%
  mutate(GOI = as.character(GOI))

data <- bind_rows(data1, data2) %>%
  mutate(GOI = parse_factor(GOI, levels = c("miR34a asRNA", "miR34a")))

#calculate ddct
groups[groups == "gene"] <- "GOI"
logical <- tibble(
  `Biological Replicate` = pull(data, `Biological Replicate`) == 2,
  `Genetic mod` = pull(data, `Genetic mod`) == "mock"
)

data <- data %>%
  ddct(., logical, groups) %>%
  folds(.)

#calculate stats
stats <- calcStats(data, `Genetic mod`, "mock", groups) %>%
  pFormat(.)

#format for plotting
lev <- c("PC3\nmock", "PC3\nmiR34a asRNA", "Skov3\nmock", "Skov3\nmiR34a asRNA", "Saos2\nmock",   "Saos2\nmiR34a asRNA")

data <- data %>%
  mutate(Condition = paste(`Cell line`, `Genetic mod`, sep = "\n")) %>%
  mutate(Condition = parse_factor(Condition, levels = lev))

stats <- stats %>%
  mutate(Condition = paste(`Cell line`, `Genetic mod`, sep = "\n")) %>%
  mutate(Condition = parse_factor(Condition, levels = lev))
#plot for rmarkdown
ggplot(data = NULL)+
  geom_violin(
    data = data,
    aes(x = Condition, y = log2fold, group = interaction(Condition, GOI)),
    trim=TRUE,
    alpha=0.5,
    fill="grey",
    colour=NA,
    scale="width",
    show.legend = FALSE
  ) +
  geom_point(
    data = data,
    aes(x = Condition, y = log2fold, colour = GOI, group = interaction(Condition, GOI)),
    size = 2, 
    alpha = 0.75,
    position = position_dodge(width = 0.9)
  ) +
  facet_grid(
    . ~ `Cell line`, 
    scales = "free", 
    space = "free"
  ) +
  geom_linerange(
    data = stats,
    aes(x = Condition, ymin = CI95l, ymax = CI95h, group = interaction(Condition, GOI)),
    position = position_dodge(width = 0.9),
    colour = "grey25",
    show.legend = FALSE
  ) +
  geom_point(
    data = stats,
    aes(x = Condition, y=mean, group = interaction(Condition, GOI)),
    colour = "black",
    size = 1,
    position = position_dodge(width = 0.9),
    show.legend = FALSE
  ) +
  geom_label(
    data = stats,
    aes(x = Condition, y = CI95h + 0.5, label = pFormat, group = interaction(Condition, GOI)),
    fill = "white",
    label.size = 0,
    label.padding = unit(0.01, "lines"),
    position = position_dodge(width = 0.9),
    show.legend = FALSE,
    fontface = "bold",
    family = "Arial Unicode MS",
    size = 3
  ) +
  scale_y_continuous("\u0394\u0394Ct")+
  labs(
    x = "Cell line",
    title = "miR34a asRNA over-expressing cell lines",
    caption = 
            "miR34a asRNA and miR34a levels in PC3, Skov3, and Saos2 cell lines stably over-expressing miR34a asRNA compared to mock over-expressing control.\n\u0394\u0394ct = House keeping gene - Target gene - mock control. Error bars indicate 95% CI and black points represent the mean."
  )+
  theme_few() +
  scale_colour_ptol() +
  theme(
    legend.position = "top",
    legend.title = element_text(size = 17),
    legend.text = element_text(size = 15),
    plot.title = element_text(
      hjust = 0.5,
      face = "bold",
      size = 20
    ),
    plot.caption = element_text(
      hjust = 0, 
      margin = margin(t = 15),
      family = "Arial Unicode MS"
    ),
    axis.title = element_text(size = 17),
    axis.title.x = element_text(margin = margin(t = 10)),
    axis.text = element_text(size = 12),
    strip.text.x = element_text(size = 17)
  )+
  guides(
    color = guide_legend(title = "Gene", override.aes = list(size = 4)))

Conclusion

miR34a asRNA over-expression causes a concominant up-regulation of miR34a in all 3 cell lines. This indicates that miR34a asRNA concordantly regulates miR34a.

Figure 3b

Introduction

miR34a is known to regulate cell cycle and, specifically, induce G1-phase arrest. Due to the increase observed in miR34a expression levels in cell lines over-expressing miR34a asRNA, we investigated if the cell cycle was perturbed in these cells.

Methods

All cell lines were cultured at 5% CO2 and 37° C with Skov3 and PC3 cells in RPMI (Hyclone) and 2 mM L-glutamine. All growth mediums were supplemented with 10% heat-inactivated FBS and 50 μg/ml of streptomycin and 50 μg/ml of penicillin. 1x10^5 cells per well in a 6-well plate and harvested after 24hrs. Cells were washed in PBS and fixed in 4% PFA at room temperature overnight. PFA was removed, and cells were resuspended in 95% EtOH. The samples were then rehydrated in distilled water, stained with DAPI and analyzed by flow cytometry on a LSRII (BD Biosciences) machine. Resulting data was further processed using Flowjo software.

Results

url <- "stableLineCellCycle.rda"
repmis::source_data(file.path(dataUrl, url))
## [1] "."
data <- .

#convert the value (precent) to log2(fraction)
data <- mutate(data, value = log2(value / 100))

#qqplot shows that the data is reasonably normal. The shapiro test is also used below. Note that a p < 0.05 indicates non-normality.
#qqnorm(data$value)

#calculate stats
stats <- data %>%
  full_join(
    filter(data, condition == "mock"), 
    by = c("Biological Replicate", "Cell line", "phase"), 
    suffix = c("", ".y")
  ) %>%
  group_by(`Cell line`, condition, phase) %>%
  summarize(
    n = n(),
    mean = (2^mean(.data$value)),
    max = (2^max(.data$value)),
    CI95l = (2^t.test(.data$value)$conf.int[1]),
    CI95h = (2^t.test(.data$value)$conf.int[2]),
    pValue = t.test(.data$value, .data$value.y)$p.value,
    shapiroTest = shapiro.test(value)$p.value
  ) %>%
  ungroup() %>%
  pFormat(.)
  #mutate(pFormat = case_when(
  #  pValue >= 0.05 & pValue != 1 ~ "n.s.",
  #  TRUE ~ pFormat
  #))

#prepare for plotting
#change value back to percent
data <- data %>%
  mutate(value = 2^value)

#calculate start and end for horizontal plot lines
stats <- stats %>%
  group_by(phase, `Cell line`) %>%
  mutate(y = max(c(.data$CI95h, .data$max))) %>%
  ungroup() %>%
  arrange(`Cell line`, phase) %>%
  mutate(
    x = rep(sort(rep(1:3, 2)), 2) - 0.25,
    xend = rep(sort(rep(1:3, 2)), 2) + 0.25
  )
#plot for rmarkdown
ggplot(data = NULL) +
  geom_violin(
    data = data,
    aes(x = phase, y = value, fill = condition),
    scale = "width",
    alpha = 0.3, 
    weight = 10,
    trim = TRUE,
    color = NA,
    show.legend = FALSE
  ) +
  geom_point(
    data = data,
    aes(x = phase, y = value, group = interaction(phase, condition), color = condition),
    position = position_dodge(width = 0.9),
    size = 2
  ) +
  facet_grid(
    . ~ `Cell line`, 
    scale = "free"
  ) +
  geom_linerange(
    data = stats,
    aes(x = phase, ymin = CI95l, ymax = CI95h, group = interaction(phase, condition)),
    position = position_dodge(width = 0.9),
    colour = "grey25",
    show.legend = FALSE
  ) +
  geom_point(
    data = stats,
    aes(x = phase, y = mean, group = condition), 
    colour = "black",
    position = position_dodge(width = 0.9),
    show.legend = FALSE
  ) +
  scale_y_continuous(labels = scales::percent_format())+
  geom_segment(
    data = stats,
    aes(x = x, y = y + 0.01, xend = xend, yend = y + 0.01),
    colour = "grey43",
    show.legend = FALSE
  ) +
  geom_label(
    data = stats,
    aes(x = phase, y = y + 0.02, label = pFormat),
    label.size = 0,
    label.padding = unit(0.01, "lines"),
    show.legend = FALSE,
    fill = "white",
    size = 3,
    family = "Arial Unicode MS"
  ) +
  labs(
    x = "Cell cycle phase",
    title = "Cell cycle",
    caption = 
            "Cell cycle analysis in miR34a asRNA over-expressing PC3, Skov3, and Saos2 stable cell lines compared to their respective mock control.\nError bars (black vertical) indicate 95% CI and horizontal lines (black) represent the mean."
  )+
  theme_few() +
  scale_fill_ptol() +
  scale_colour_ptol(labels = c("PC3 mock", "PC3 miR34a asRNA")) +
  theme(
    plot.title = element_text(
      hjust = 0.5,
      face = "bold",
      size = 22
    ),
    plot.caption = element_text(
      hjust = 0, 
      margin = margin(t = 15),
      family = "Arial Unicode MS"
    ),
    legend.position = "top",
    axis.title = element_text(size = 17),
    axis.title.x = element_text(margin = margin(t = 10)),
    axis.text = element_text(size = 15),
    strip.text.x = element_text(size = 15),
    axis.title.y = element_blank(),
    legend.title = element_text(size = 17),
    legend.text = element_text(size = 15)
  )+ 
  guides(color = guide_legend(title = "Cell line")) 

Conclusion

An increase in the percentage of cells in G1-phase can be seen in cell lines overexpressing miR34a asRNA. In addition, a concominant decrease of cells in G2-phase can be observed.

Figure 3c

Introduction

As p53 is a well-known regulator of senescence, growth, and apoptosis and has been shown to specifically important under conditions of cellular stress, such as starvation. Therefore, we investigated whether overexpression of miR34a AS affects growth rate under normal and starvation conditions by measuring changes in confluency over time.

Methods

10^4 PC3 cells, either miR34a asRNA overexpressing or WT, were seeded in 96 well plates. After attachment (3-4h) media was replaced with either RPMI (Gibco, life technology) (supplemented with 2mM L-glutamine, 50ug/ml Penicillin-Streptomycin and 10% Fetal Calf Serum) or HBSS. Cells were incubated in Spark Multimode Microplate reader for 48 hours at 37°C with 5% CO2 in a humidity chamber. Confluency was measured every hour.

Analysis

Load the data.

url <- "growthStarvation.rda"
repmis::source_data(file.path(dataUrl, url))
## [1] "."
data <- .

#normalize to first time point for control (Mock)
.normFun <- function(
  data,
  Confluency,
  Treatment,
  Biological.Replicate,
  Technical.Replicate,
  Cell.line
){
  bool1 <- data$Treatment == Treatment
  bool2 <- data$`Biological Replicate` == Biological.Replicate
  bool3 <- data$`Technical Replicate` == Technical.Replicate
  bool4 <- data$Time == 0
  bool5 <- data$`Cell line` == Cell.line
  bool <- bool1 & bool2 & bool3 & bool4 & bool5
  Confluency / pull(data, Confluency)[bool]
}

data <- data %>%
  group_by(Treatment, `Biological Replicate`, `Technical Replicate`, Time, `Cell line`) %>%
  mutate(
    normConfluency = .normFun(
      data,
      Confluency,
      Treatment,
      `Biological Replicate`,
      `Technical Replicate`,
      `Cell line`
    )
  ) %>%
  ungroup()

#calculate mean for each technical replicate
data <- data %>%
  group_by(Treatment, `Biological Replicate`, Time, `Cell line`) %>%
  summarize(
    techMean = mean(normConfluency)
  ) %>%
  ungroup()

#calculate stats
stats <- data %>%
  filter(Time != 0) %>%
  group_by(`Cell line`, Treatment, Time) %>%
  summarize(
    min = min(techMean),
    max = max(techMean),
    mean = mean(techMean),
    median = median(techMean),
    CI95l = t.test(techMean)$conf.int[1],
    CI95h = t.test(techMean)$conf.int[2]
  )

Results

#plot for Rmarkdown
ggplot(data = NULL) + 
  geom_ribbon(
    data = stats,
    aes(x = Time, ymin = CI95l, ymax = CI95h, group = `Cell line`, fill = `Cell line`), 
    alpha = 0.2,
    show.legend = TRUE
  ) +
  geom_linerange(
    data = stats,
    aes(x = Time, colour = `Cell line`, ymin = min, ymax = max), 
    position = position_dodge(width = 1),
    show.legend = FALSE
  ) + 
  geom_point(
    data = stats,
    aes(x = Time, y = median, group = `Cell line`, colour = `Cell line`), 
    position = position_dodge(width = 1),
    show.legend = FALSE
  ) +
  facet_grid(. ~ Treatment) +
  scale_x_discrete(breaks = seq(0, 47, 5)) +
  theme_few() +
  scale_fill_ptol() +
  scale_colour_ptol() +
  theme(
    plot.title = element_text(
      hjust = 0.5,
      face = "bold",
      size = 20
    ),
    plot.caption = element_text(
      hjust = 0, 
      margin = margin(t = 15),
      family = "Arial Unicode MS"
    ),
    legend.position = "top",
    legend.title = element_text(size = 17),
    legend.text = element_text(size = 15),
    axis.title = element_text(size = 17),
    axis.title.x = element_text(margin = margin(t = 10)),
    axis.text = element_text(size = 15),
    strip.text.x = element_text(size = 17)
  ) +
  labs(
    x = "Time (hours)", 
    y = "Normalized % confluency",
    title = "miR34a asRNA's role in growth regulation.",
    caption = "The effects of miR34a asRNA overexpression on growth in normal and starvation conditions in the PC3 prostate cancer cell line."
  ) +
  guides(
    fill = guide_legend(title = "Cell line"), 
    colour = guide_legend(title = "Cell line")
  )

Conclusions

miR34a asRNA overexpression causes a minor decrease in cell growth under normal conditions (RPMI) although under conditions of cellular stress through starvation, this effect is increased dramatically.

Figure 3d

Introduction

miR34a asRNA over-expression was shown to increase the expression levels of miR34a. We investigated if this increase of miR34a was due to increased levels of transcription at the miR34a promoter.

Methods

All cell lines were cultured at 5% CO2 and 37° C with Skov3 and PC3 cells in RPMI (Hyclone) and 2 mM L-glutamine and Saos2 cells cultured in DMEM high modified (Hyclone). All growth mediums were supplemented with 10% heat-inactivated FBS and 50 μg/ml of streptomycin and 50 μg/ml of penicillin. Cells were plated at 2-3x10^6 cells in a 15cm dish and incubated overnight. Cells were cross-linked in 1% formaldehyde, quenched in glycine (0.125M), resuspended in cell lysis buffer (5mM PIPES, 85mM KCL, 0.5% NP40, protease inhibitor), and sonicated (Bioruptor Sonicator Diagenode) in nucleus lysis buffer (50mM TRIS-HCL pH 8.0, 10mM EDTA, 1% SDS, protease inhibitor). The samples were incubated overnight at 4 °C with the appropriate antibody. Salmon sperm DNA/Protein A–agarose (Upstate/Millipore) was used to pull down the antibody. DNA was eluted in elution buffer (1% SDS, 100 mM NaHCO3), followed by reverse cross-linking, RNase-A and protease-K treatment. The DNA was eluted by using a Qiagen PCR purification kit. The following antibodies were used (4 μg/sample): phospho-polymerase II (abcam: ab5095). QPCR was carried out using KAPA 2G SYBRGreen (Kapa Biosystems) using the Applied Biosystems 7900HT machine with the cycling conditions: 95 °C for 3 min, 95 °C for 3 s, 60 °C for 30 s.

Primers
primers <- data.frame(
    name=c(
        "miR34aChIP F1",
        "miR34aChIP R1"
    ),
    sequence=c(
       "AAA GTT TGC AAA GAA GGA GGC GGG",
       "AGG GAA GAA AGA ACT AGC CGA GCA"
    )
)

Results

url <- "stableLinePolIIChIP.rda"
repmis::source_data(file.path(dataUrl, url))
## Downloading data from: https://github.com/GranderLab/miR34a_asRNA_project/raw/master/data/stableLinePolIIChIP.rda
## SHA-1 hash of the downloaded data file is:
## f8fff1b9860dd20d29ad0a08060b4cf18afd0938
## [1] "."
data <- .

#calculate means of technical replicates
data <- mutate(data, techMean = rowMeans(select(data, Quantity1, Quantity2)))

#calculate percent of input
.percentOfInput <- function(data, techMean, Biological.Replicate, condition) {
  bool1 <- data$`Biological Replicate` == Biological.Replicate
  bool2 <- data$condition == condition
  bool3 <- data$sample == "input"
  bool <- bool1 & bool2 & bool3
  techMean / pull(data, techMean)[bool]
}

data <- data %>%
  group_by(`Biological Replicate`, condition, sample) %>%
  mutate(
    percentInput = .percentOfInput(data, techMean, `Biological Replicate`, condition),
    logPercentInput = log2(percentInput)
  ) %>%
  ungroup() %>%
  filter(sample != "input")

#calculate fold over control (mock)
.foldOverControl <- function(data, logPercentInput, Biological.Replicate) {
  bool1 <- data$`Biological Replicate` == Biological.Replicate
  bool2 <- data$condition == "mock"
  bool <- bool1 & bool2
  logPercentInput - pull(data, logPercentInput)[bool]
}

data <- data %>%
  group_by(`Biological Replicate`) %>%
  mutate(
    logfoldOverControl = .foldOverControl(data, logPercentInput, `Biological Replicate`),
    foldOverControl = 2^logfoldOverControl
  )

#calculate stats 
stats <- data %>%
  group_by(condition) %>%
  summarize(
    n = n(),
    mean = mean(logfoldOverControl),
    CI95l =t.test(logfoldOverControl)$conf.int[1],
    CI95h =t.test(logfoldOverControl)$conf.int[2],
    pValue = t.test(logfoldOverControl, mu = 0)$p.value
  ) %>%
  pFormat(.)
#plot for Rmarkdown
ggplot(data = NULL)+
  geom_violin(
    data = data,
    aes(x = condition, y = logfoldOverControl, fill = condition),
    trim = TRUE,
    alpha = 0.3,
    colour = NA,
    show.legend = FALSE
  ) +
  geom_point( 
    data = data,
    aes(x = condition, y = logfoldOverControl, colour = condition),
    size = 4
  ) +
  geom_linerange(
    data = stats,
    aes(x = condition, ymin = CI95l, ymax = CI95h),
    colour = "black",
    show.legend = FALSE
  ) +
  geom_point(
    data = stats,
    aes(x = condition, y = mean),
    colour = "black",
    show.legend = FALSE
  ) +
  geom_segment(
    data = stats,
    aes(x = 1, y = CI95h + 0.25, xend = 2, yend = CI95h + 0.25),
    colour = "grey43"
  ) +
  geom_label(
    data = stats,
    aes(x = 1.5, y = CI95h + 0.35, label = pFormat),
    label.size = 0,
    label.padding = unit(0.01, "lines"),
    show.legend = FALSE,
    fill = "white",
    size = 5,
    family = "Arial Unicode MS"
  ) +
  scale_y_continuous("log2(Fold)") +
  labs(
    x = "Cell line",
    y = "log2(Fold polII enrichment)",
    title = "Phosphorylated Poll II ChIP",
    caption = 
            "Phosphorylated Poll II ChIP in miR34a asRNA over-expressing PC3 stable cell lines compared to their respective mock control.\nError bars (black vertical) indicate 95% CI and black points represent the mean."
  )+
  theme_few() +
  scale_fill_ptol() +
  scale_color_ptol() +
  theme(
    plot.title = element_text(
      hjust = 0.5,
      face = "bold",
      size = 20
    ),
    plot.caption = element_text(
      hjust = 0, 
      margin = margin(t = 15),
      family = "Arial Unicode MS"
    ),
    legend.position = "none",
    axis.title = element_text(size = 17),
    axis.title.x = element_text(margin = margin(t = 10)),
    axis.text = element_text(size = 15)
  )

Conclusion

Phosphorylated pol-II was increased at the miR34a promoter in miR34a asRNA over-expressing cell lines indicating that the regulation of miR34a by miR34a asRNA is transcriptional.

Session info

sessionInfo()
## R version 3.4.2 (2017-09-28)
## Platform: x86_64-apple-darwin15.6.0 (64-bit)
## Running under: macOS High Sierra 10.13
## 
## Matrix products: default
## BLAS: /Library/Frameworks/R.framework/Versions/3.4/Resources/lib/libRblas.0.dylib
## LAPACK: /Library/Frameworks/R.framework/Versions/3.4/Resources/lib/libRlapack.dylib
## 
## locale:
## [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
## 
## attached base packages:
## [1] stats     graphics  grDevices utils     datasets  methods   base     
## 
## other attached packages:
##  [1] bindrcpp_0.2                  miR34AasRNAproject_0.0.0.9000
##  [3] RCurl_1.95-4.8                bitops_1.0-6                 
##  [5] ggthemes_3.4.0                printr_0.1                   
##  [7] dplyr_0.7.4                   purrr_0.2.4                  
##  [9] readr_1.1.1.9000              tidyr_0.7.2                  
## [11] tibble_1.3.4                  ggplot2_2.2.1                
## [13] tidyverse_1.1.1              
## 
## loaded via a namespace (and not attached):
##  [1] tidyselect_0.2.2    reshape2_1.4.2      haven_1.1.0        
##  [4] lattice_0.20-35     colorspace_1.3-2    htmltools_0.3.6    
##  [7] yaml_2.1.14         rlang_0.1.2         R.oo_1.21.0        
## [10] R.utils_2.5.0       foreign_0.8-69      glue_1.1.1         
## [13] modelr_0.1.1        readxl_1.0.0        R.cache_0.12.0     
## [16] bindr_0.1           plyr_1.8.4          stringr_1.2.0      
## [19] munsell_0.4.3       gtable_0.2.0        cellranger_1.1.0   
## [22] R.methodsS3_1.7.1   rvest_0.3.2         psych_1.7.8        
## [25] evaluate_0.10.1     labeling_0.3        knitr_1.17         
## [28] forcats_0.2.0       curl_3.0            parallel_3.4.2     
## [31] repmis_0.5          highr_0.6           broom_0.4.2        
## [34] Rcpp_0.12.13        scales_0.5.0        backports_1.1.1    
## [37] jsonlite_1.5        mnormt_1.5-5        hms_0.3            
## [40] digest_0.6.12       stringi_1.1.5       grid_3.4.2         
## [43] rprojroot_1.2       tools_3.4.2         magrittr_1.5       
## [46] lazyeval_0.2.0      pkgconfig_2.0.1     data.table_1.10.4-2
## [49] xml2_1.1.1          lubridate_1.6.0     assertthat_0.2.0   
## [52] rmarkdown_1.6       httr_1.3.1          R6_2.2.2           
## [55] nlme_3.1-131        compiler_3.4.2